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This routine is a subroutine version of M 2k » It solves 

the matrix equation AX = B given A and B where 

1. A is a non-singular square matrix of order n, 

2. B is of size (n x m). B may be the identity matrix 
I of size (n x m). 

The magnitudes of n and m are governed by the storage 
available in the Williams memory and on the drum. The 
matrix uses 

2n + km + 1 
storage locations in the Williams memory and a minimum of 

(B = I) 

or n(n + m + 2) , (B f i) 



n(n +l) , / , , \ 
-i-75 <- + n(m + 1) , 



storage locations on the drum. 

The subroutine takes the matrices (or matrix) from the 

drum and puts the solution back on the drum. 

The parameter S3 is used to determine the first location, 

S, of Williams memory available for storage. It then 

uses the 2n + km + 1 successive locations starting with 

S. The parameter has the form 

00 F 

00 SF 
This subroutine assumes the equations to be solved are 
on the drum, beginning at location d Q , stored by row 
with a sum check after each row, in the form 



a ii' a i2' a iy 

b il' b 12' " 
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SUM CHECK 

SUM CHECK 

* 
• 

SUM CHECK 
SUM CHECK 
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where a . and b. . are the elements of the matrices A 
and B. In the case of inversion B is not placed on the 
drum and the successive rows of A follow one another 
immediately. 

The routine uses ~p *~ + n(m + l) drum locations for 

temporary storage starting with location d 1 ., 

The solution is stored by row with a sum check after 

each row starting with drum location d . The scaling 

factors are as computed in M 2k and are stored as an 

extra row of the solution matrix. This takes (n + l)(m + l) 

locations. 

Before entry to this subroutine, Williams memory locations 

3 through 7 must contain 

3: n, order of matrix A 

4: m, number of columns in matrix B (m ■ n 
for inversion, B = I) 

5: d n , as defined above 



6: d^ " 

7:*,, " 
Entry is of the form 


n n 


p 50 F 




50 pF 
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for the solution of the 


linear matrix equation, and 


p JO F 




50 P F 





26 (M27) 
for inversion. Control is returned to the right side 
of location p + 1. 

If a solution has been obtained, location contains 
zero. If A is singular, hence no solution, location 
is set to l/2. 
DURATION: Characteristic times are 

n, m = 6 7 seconds 

n, m = 10 21 seconds 

n, m = 15 .5^ seconds 
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NOTE! 



For large n, the time goes up as n . Due to the large 
use of the drum the routine should only he used for 
n > 20. M Ik may he used for lower order matrices. 

1. Subroutine Y 1 is contained in the subroutine itself 
and may he used separately. It has the symbolic address 

Of i). 

2. In the case of inversion, care should he taken to keep 
d Q and d, far enough apart. One must have 

dj_ < d Q - [n(m + 1) - S^Lii ] 

or d 1 > d Q + n(n + l). 

Suitable- choices of d , d, , and d will allow one to 
have the original matrices availahle after obtaining a 
solution. 

3. A stop with 

FF 010 

k-2 

signifies a sum check stop in Y 1. A white switch 
start will attempt to read again. 
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1 


1*2 204L 
40 8F 




1 Presets 1 
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L5 5F 




I 
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L4 21% 








40 29L 




d o 


3 


L5 6f 
l4 215L 






4 


40 183L 
L5 7F 




d ! 


5 


L4 215L 








40 190L 




d 2 


6 


L5 3F 
00 20F 






7 


1^6 30L 




n 




L5 4f 






8 


50 205L 
00 20F 
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46 42L 
46 191L 




m 


10 


46 202L 
L4 30L 






11 


46 210L 
46 184L 




M= n + m 
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L5 210L 
L4 3F 






13 


42 122L 
42 143L 




y = S3 + n 
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00 20F 
46 40L 
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46 122L 
46 13TL 






16 


46 176L 








46 189L 
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L5 122L 
L4 4f 






18 


42 60L 
42 65L 




t=y + n»S3+/f 


19 


42 66L 
42 211L 






20 


42 214L 
00 20F 
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46 56L 
46 95L 






22 


46 131L 
46 162L 






23 


46 212L 
L4 210L 
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46 188L 
L4 42L 




u = S3 + 2fi 


25 


46 200L 
LO 42L 




v « S3 + 3/ 1 


26 


10 20F 




u , i 




42 213L 
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28 


41 6F 
4l 7F 
50 S3 
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50 28L 
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k Call row r of A from Drum 


30 


00 (n)F 
F5 29L 
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L4 3F 
40 29L 






32 


4o'4lL 






33 


L5 6T 
36 40L 
L5 122L 




»B / It 


> 


42 35L 
L4 6F 
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HOPES PACK 5 


55 


k2 58L 
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kl (y')F 


* 


Augment A with B = I 


56 


F5 55L 
42 55L 


r 




57 


LO 2l4L 
5255L 






58 


L5 209L 








kO (y + r)F 


35 




59 


26 44L 








00 F 




Waste 


-■■ko 


50 (y)F 
50 40L 


14* 




kl 


26 (Yl) 
00 (d)F 


r 


Call rov r of B from Drum 


42 


00 (m)F 
F5 4lL 






k3 


L4 kF 
kO 29L 






44 


L5 185L 




m 




M 


kO 57L 

ko 96L 

L5 210L 
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Preset Drum 
" Compaands 


k6 


46 58L 
k6 9TL 








kf 


L5 210L 








L4 7F 




x = S5 + i 


48 


42 64L 
42 69L 






k9 


42 75L 
42 89L 






50 


42 91L 
42 114L 




-Set x addresses 


51 


00 20F 
46 61L 






52 


46 65L 




, 




46 67L 
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66 
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46 75L 




46.75L 




46 U3L 




L3 6F 




&6 7P 




36106L 




50 (t)F 


21 


50 56L 




26 (Yl) 
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00 (d)F 
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00 (v)F 


46 


26 59L 




L5 206L 




46 70L 




46 71L 




L3 (t)F 


18 


l£ (x)F 


51* 


32 64L 




47 70L 




50 205L 




L5 (x)F 


52 


66 (t)F 


18« 


26 68L 




50 (x)F 


48 


S3 F 




32 67L 




50 208L 




75 (t)p 


19 


66 (x)f 
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47 71L 




4l 5F 




SI F 




40 2F 




L5 (x)F 


48- 


40 (l^f 
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L5 (t)F 
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80 

62 
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-i > ^ 



Call row i for 
eliioinatioji 



-Beset interchange 



x > tt 



1^ ->0 



*k = 



¥ 



-x = 0? 



X 



1^ ->o 
Clear M "box 
Store - k 



-Interchange 



LOCATION 


ORDER 




HOIKS 


PACE 5 


71 


ho (I^)F 
50 2F 


6067* 


j 




72 

73 


7J IF 
L4 F 
40 (x)F 
L3 (x)F 


3 "' , 


i ~ t. 1 

t, - kx,, k - — 


or 


7^ 


L6 5F 


■ 


-M>|x J M 
J 






56 76L 
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L7 (x)F 
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"-Replace |M| by |x ' | 

J 






40 5F 
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L5 IF 








77 


kO (t)F 
F5 76L 
k2 76L 


77* 103 


Complete change 




78 


1*2 70L 
L5 73L 




Advance 




79 


IA 206L 
40 75L 




addresses 




80 


k2 69L 
1*6 75L 




/ 




81 


LO 212L 
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52 69L 
L5 207L 
kO 90L 




More elements! 




83 


LL 5F 








52 85L 






Qk 


L5 69L 
1*6 90L 
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26 89L 
L5 90L 1 
















86 


k-2 90L 










L5 5F 
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00 IF 
1*0 5F 








88 


LL 5F 




Scale to 
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32 85L 
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91 
92 

95 
9^ 
95 
96 

97 
98 

99 
100 
101 
102 

103 
104 

105 
106 



50 205L 
L5 (x)F 
10 (l)F 
00 (l)F 
50 205L 
40 (x)F 
P5 91L 
42 91L 
42 89L 
LO 211L 
52 89L 
26 95L 
JO (t)F 
50 95L 
26 (Yl) 
00 (d)F 
00 (v)F 
L5 97L 
10 20F 
F4 96L 
40 96L 
40 57L 
L5 97L 
LO 206L 
46 97L 
46 58L 
L5 211L 
42 70L 
42 76L 
L5 212L 
46 95L 
F5 7F 
42 7F 
26 47L 
L3 7F 
L6 6f 
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8k 
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95 
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110 f 
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99 

46' 101 
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-Put row i back on dtfum 



Reset t address of 



i ->i + 1 
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109 
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111 

112 

113 

114 

115 
116 
117 
118 

119 
120 
121 
1A 

124 



ORDER 

36 110L 
F5 6F 
42 6f 
LO 3F 
36 121L 
22 27L 
L5 113L 
46 95L 
L3 6F 
32 112L 
26 95L 

41 5F 
L3 (x)F 
L6 5F 
32 115L 
L7 (x)F 
40 5F 
F5 114L 

42 114L 
00 20F 
46 113L 
F5 7F 
42 7F 
LO 3F 
LO 4F 
36 82L 
26 113L 
00 F 

4l 5F 
L5 209L 
40 (y)F 
L3 (y)F 
32 203L 
L5 96L 
40 132L 
L597L 



»OTBB J 
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r -»r + 1 
r - «T 



117 54 



116 • 50 « 



Prepare to put row r 
back on Drum 



15 • 
13 



Clear k 

Stores 
►•Sing? S < 2 



Back Substitution , 

i — 
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Set Drum Addresses 
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127 
128 

150 
131 
132 
133 
134 
135 
136 

137 
138 

139 
140 
l4i 

142 



ORDER 

46. 133L 
4l 6T 
kl 7F 
J& 132L 
JO 4f 
FO 6F 
40 132L 
40 165L 
L5 3^3L 
L4 206L 
46 133L 
46 164L 
50 (t)F 
50 131L 
26 (Yl) 
00 (d)F 
00 (v)F 
L5 211L 
F4 6¥ 
42 l40L 
L4 5F 
42 137L 
L5 137L 
46 140L 
50 (y)F 

71 (b)F 
40 F 

22 144L , 
00 F 
S5 F 
50 (x)F 
74 (a)F 
L4. F 
40 F 
LL F 
32 144L 






"B^T 



22 



)l24 

J 128 

125 130 



15 f 
135 



136 

134 



:} 



145' 
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Clear Counters 



Adjust Druft Addresses 



—Call row (n - r) froa 



Reset addresses 



Compute 

- n 

z 

r=Ul 



a J x 

lr r 
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1*3 


50 2091. 




7J (y)y 


144 


26 122k 




15 140L 


1*3 


L0 206L 




40 140L 


146 


42 150L 




42 155L 


147 


46 155L 




F5 7F 


148 


42 7F 




L5 6F 


149 


10 7? 




32 139L 


150 


4l 7F 




L3 (a^jF 


151 


32 203t 




l£ F 
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36BT3L 




26 153L 


155 


L5 F 




66 (a i:L )F 


154 


22 154L 
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40 (*)F 




L5 5F 
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32 164L 




L5 155L 


157 


L4 210L 




L4 210L 
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46 161L 




L5 211L 
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L4 5F 




B4 6F 
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42 161L 




26 16IL 
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Rescale 
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Advance addresses 
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for rov r 



Reset 1 
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DiYiaion bad? 
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1 



k = 0? 



Restore row (n - r) 
of augmented matrix 
to Drum vlth 

Vl, r 
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40 (o<)F 
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30 162L 
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I63 
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00 (d)F 
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164 


00 (y)F 
F5 6F 


130' 


, 


165 


42 6f 
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LO 3F 




Count n rows 
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36 I67L 








22 126L 
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10 213L ' 
40 169L 






168 


4l F 

26 I69L 
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L5 S3 


1 167' 






ho (u)F 


J171 




170 


L5 169L 
L4 206L 






171 


4o 169L 

F5 F 




Store eoliun 
■" k for later use 


172 


42 F 
LO 4F 






173 


36 174L 

26 169L 
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L5 200L 
10 20F 






175 


L4 5F 
42 176L 






176 


L5 (3T)F 


16 

1T5* = 




177 


40 (v 1 )F 
F5 5F 


1 




42 5F 




k -*k + 1 


178 


LD 4f 
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32 179L 
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22 121L 




J 
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L5 211L 
W 206L 
42 188L 
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4l 2F 
26 182L 
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50 S3 








50 182L 
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26 (Yl) 
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00 (v)F 
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L5 184L 
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10 20F 
F4 183L 
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4o 183L 








10 184L 
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LO 206L 
46 184L 




Arrange 
by ro* 
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L5 (u)F 


24 196 






4o (v)f 


180 ' 197 




189 


JO (y)F 
50 I89L 


l6» 192' 




190 


26 (Yl) 
00 (d)F 


J 194 




191 


00 (m)F 
L5 189L 
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LO 206L 
46 189L 
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L5 190L 

f4 4f 
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40 190L 
40 201L 
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L5 188L 
L4 206L 
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46 188L 








LO 205L 
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wtm '»arnr 
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kz 188L 
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F5 2F 
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k2 2F 








LO 5P 
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36 200L 








26 182L' 
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JO (y)F 
50 200L 


25 
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26 (Yl) 
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00 (d)F 


J 19*' 


row n + 1 of x 


202 


00 (m)F 


10 






Hi F 
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22 2C*L 








k9 F 




s- If singular 


20V 


22 204L 
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o« 


Exit 
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00 F 
00 2F 
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00 IF 
00 IF 
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10 IF 
00 IF 
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7L H095F 

LL H095F 
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00 F 
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00 (/<)F 

00 S3 


n 




211 


JO 205L 








40 (t)F 


19 r 




212 


MO (t)F 

L3 F 
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213 


L5 S3 








40 (u)F 
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12 38L 
41 (t)? 
26 (II) 
00 F 
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